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art  softvare  has  been  created  for  the  nonparametric  estima¬ 
tion  and  contour  display  of  densities  of  dimensions  through  five.  A  non¬ 
parametric  mode  finding  algorithm  has  been  developed  for  densities  of 
dimension  up  to  one  hundred.  An  algorithm  for  data  based  simulation  from 
high  dimensional  distributions  has  been  created  and  employed  on  behind 
armour  data.  In  the  related  area  of  parameter  estimation  of  stochastic 
processes,  simulation  based  algorithms  have  been  developed  which  enable 
the  user  to  proceed  directly  from_the  axioms  of  the  process  nlus  data  to 
the  estimation  of  parameters.  T  ■ 
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Statement  of  Problem  Studied 


I 


Our  task  has  been  to  develop  algorithms*  and*  to  whatever  extent 
possible*  portable  software  for  the  data  based  nonparametric  estimation 
of  densities  of  high  dimensionality.  Ve  hoped  to  employ  graphical  and 
interactive  techniques  whenever  possible.  In  cases  where  the  dimen¬ 
sionality  was  very  high*  we  hoped  to  develop  an  algorithm  for  the  find¬ 
ing  of  modes  which  might  then  be  used  as  origins  for  further  investiga¬ 
tions.  In  addition*  we  hoped  to  examine  related  problem  areas*  such  as 
parameter  estimation  for  applied  stochastic  process  models.  It  was  hoped 
that  interaction  with  the  Ballistics  Research  Laboratory  at  Aberdeen 
would  be  fruitful.  Also*  it  was  hoped  that  the  Texas  Medical  Center  and 
NASA  would  provide  useful  local  data  sources  on  which  to  test  our  algo¬ 
rithms  . 
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Summary  of  Moat  Important  Results 


The  main  thrust  of  our  vork  has  been  in  the  nonparametric  estima¬ 
tion  and  representation  of  densities  of  dimensions  of  three  and  greater. 
An  algorithm  for  multicolour*  time  lapsed  density  contour  display  has 
been  developed  by  Scott  and  is  nov  portable.  A  16mm  film  ha6  been  made 
of  the  application  of  this  technique  to  estimate  and  display  densities 
of  dimensions  three*  four  and  five.  This  film  has  been  shown  at  the 
Aberdeen  Proving  Ground*  at  the  1983  DOD  Conference  at  the  USAF  Academy* 
at  the  ONR  1983  Luray  Conference,  in  addition  to  some  dozen  nonDOD 
conferences.  The  principal  data  sets  used  in  explication  of  the  Scott 
algorithm  are  ground  observations  from  satellite  overflights. 

The  nonparametric  estimation  and  representation  of  densities  of 
very  high  dimensionality  (up  to  one  hundred)  is  a  problem  addressed  in 
the  doctoral  dissertation  of  Steven  Boswell*  presently  of  Lincoln 
Laboratories.  Boswell's  algorithm  seeks  out  the  modes  of  the  density 
which  generated  a  data  set.  These  modes  are  then  to  be  used  as  origins 
for  further  analysis. 

The  dissertation  of  Hathaway  uses  the  EM  algorithm  to  address  the 
problem  of  estimating  a  density  as  a  mixture  of  Gaussian  densities. 
Many  of  the  instablity  problems  generally  encountered  with  mixture  tech¬ 
niques  are  eliminated. 

The  problem  of  using  a  data  base  as  a  vehicle  for  simulation  has 
been  addressed  by  Taylor  and  Thompson  using  behind  armour  data  from  the 
Aberdeen  Proving  Ground.  Their  procedure  gives  the  same  effects  as  those 
which  would  have  been  obtained  if  the  underlying  density  had  first  been 
estimated  nonparametrically  and  then  used  to  build  a  random  number  gen¬ 
erator.  However*  by  eliminating  the  density  estimation  step*  the 
Taylor-Thompson  algorithm  provides  enormous  savings  both  in  computer 
running  times*  and*  more  importantly*  in  the  ease  of  use. 

The  development  of  simulation  techniques  for  parameter  estimation 
in  complex  stochastic  process  situations  has  been  a  major  activity  for 
Thompson  and  his  colleagues  at  the  Texas  Medical  Center*  Atkinson*  Bar- 
toszynski  and  Brown.  Not  only  are  these  techniques  more  stable  and 
quicker  than  those  based  on  approximation  theoretic  procedures*  but* 
since  they  require  only  the  axioms  of  the  process*  rather  than  some 
"closed  form  solution*”  they  are  much  easier  to  use.  A  problem  which 
had  required  a  man  year  to  code  using  stsndard  techniques  was  coded  in 
one  week  using  the  simulation  based  strategy.  The  development  of  this 
algorithm  enables  the  use  of  stochastic  models  in  real  world  situations 
which  earlier  had  defied  such  applied  modeling. 
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